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Genomic imprinting, the preferential expression of maternal or paternal alleles of imprinted 
genes, is often maintained through expression of imprinted long non-coding (Inc) "anti- 
sense" RNAs. These may overlap imprinted transcripts, and are expressed from the 
opposite allele. Previously we have described brain region-specific imprinted expression of 
the Dio3 gene in rat, which is preferentially modified by fetal ethanol exposure. The Dio3os 
(opposite strand) transcript is transcribed in opposite orientation to Dio3 in mouse and 
human, partially overlaps the Dio3 promoter, and mirrors total Dio3 developmental expres- 
sion levels. Here, we present that the rat Dio3os transcript(s) exhibits brain region-specific 
imprinted expression patterns similar to those of Dio3. Rat Dio3os transcript expression is 
also similarly modified by fetal ethanol exposure. Uniquely, both Dio3 and Dio3os expres- 
sion occur on the same, rather than opposite, alleles, as determined by strand-specific 
RT-PCR. Future studies will require direct manipulation of the Dio3os transcript to deter- 
mine whether the novel paralleling of total and allele-specific expression patterns of this 
sense/antisense imprinted gene pair reflects an as-yet undefined regulatory mechanism 
for IncRNA mediated tissue-specific imprinted expression, or rather is a consequence of a 
more straightforward, but previously undescribed transcriptional coregulation process. 
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INTRODUCTION 

Over the past several years, long non-coding RNAs (IncRNAs) have 
been found to be transcribed across the genome (e.g., Katayama 
etal., 2005; Huang etal., 2011). Many of these transcripts have a 
role in the regulation of gene expression, such as via direct overlap 
with genes or their promoters (reviewed in Amaral and Mattick, 
2008) or interaction with and modification of epigenetic chro- 
matin marks (Khalil etal, 2009). Historically, the first IncRNAs 
were identified within regions of genomic imprinting, charac- 
terized by preferential expression of maternal or paternal alleles 
of imprinted genes. Such imprinted expression is often main- 
tained through expression of "antisense" IncRNAs, which overlap 
imprinted transcripts and are expressed from the opposite parental 
allele, potentially "blocking" expression from the "sense" allele 
on this parental chromosome (reviewed in Peters and Robson, 
2008). In addition to direct "antisense" overlap, several imprinted 
IncRNAs have also been shown to exhibit long-distance cis effects 
on other genes within the regulatory clusters through interaction 
with chromatin and DNA modifying proteins (Nagano et al., 2008; 
Pandey etal., 2008), similar to what is being shown for IncRNAs 
in non-imprinted regions. 

The type 3 deiodinase gene, Dio3, lies at the distal end of 
a large cluster of imprinted genes. Dio3 is preferentially pater- 
nally expressed in most tissues (Charalambous and Hernandez, 
2012), as are the three other imprinted coding genes in this 
cluster (Hagan etal, 2009). In contrast, but consistent with the 
emerging realization of complex allele-preferential gene expres- 
sion patterns across the genome, we have recently described brain 
region-specific imprinted expression of Dio3 in rat (Sittig etal, 



2011a). This imprint profile is preferentially modified by fetal 
ethanol exposure and correlates with behavioral alterations (Sittig 
etal., 2011b), providing the first evidence of functional conse- 
quences of brain region-specific imprinted expression profiles. 
Although the mechanisms underlying these complex imprint pat- 
terns are not yet understood, many possibilities arise through 
comparison with other imprinted loci, as described above. Imme- 
diately adjacent to the Dio3 gene is a IncRNA transcript, Dio3os, 
that is transcribed in opposite orientation to Dio3 in mouse and 
human, partially overlaps the Dio3 promoter (Hernandez, 2005), 
and mirrors total Dio3 developmental and diurnal (Labialle etal., 
2008) expression levels. Unlike the large number of large and small 
non-coding RNA genes within this cluster which are maternally 
expressed (daRocha et al, 2008), Dio3os has been reported to have 
biallelic, rather than imprinted, expression in murine embryos 
(Tierling etal., 2006) and in adult mouse cortex (Labialle etal., 
2008). To investigate whether the Dio3os transcript exists in rat 
and thus might contribute to complex Dio3 expression regulation, 
we have analyzed the total and allele-specific expression patterns 
of rDio3os in naive and ethanol-exposed animals. These analyses 
demonstrate a combination of features and complexities of rDio3 
expression that together define a novel category within imprinted, 
IncRNAs. 

RESULTS 

CHARACTERIZATION AND BRAIN REGION-SPECIFIC IMPRINTING OF 

rDio3os 

Transcripts across the upstream rDio3 region were identified by the 
single poly-A-containing expressed sequence tag (EST) contained 
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in the NCBI database, and by RT-PCR and rapid amplification 
of cDNA ends (RACE) across the region from placenta and both 
neonatal and adult rat frontal cortex (FX) and hippocampus (HP). 
Both unspliced and alternatively spliced transcripts were identi- 
fied. rDio3os exons overlap those of mouse (Tierling et al., 2006) 
and extend into the 5' rDio3 GC rich region, but splicing and 
exon-intron boundaries are not identical. Mouse and rat show 
an average of 83% identity across overlapping regions. Alterna- 
tively spliced RT-PCR products are indicated in Figure 1A, as are 
locations of representative primer pairs used for strand-specific 
RT-PCR (below). We have not been able to identify rDio3os splice 
variants that extend through and overlap the rDio3 5'UTR and 
transcription start site, as has been reported in mouse but not 



humans (Hernandez etal, 2004), although strand-specific RT- 
PCR has identified the presence of at least an 154 bp opposite 
strand (OS) transcript 60 bp upstream of the standard rDio3 
5'UTR, within the minimal promoter and potentially overlapping 
an alternative transcription start site identified in keratinocytes 
(Dentice et al., 2007; Figure IB, "Z"). 

UPSTREAM EXPRESSION IS DERIVED EXCLUSIVELY FROM THE 
OPPOSITE STRAND TO rDio3 

The EST BI274690, ~3-4 kb upstream of rDio3 (inclusive), does 
not encode an open reading frame but contains a poly-A sequence, 
strongly suggesting that its direction of transcription is opposite 
to that of rDio3 and that it represents the rDio3os transcript rather 
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FIGURE 1 | (A) Genomic organization of rDio3os transcripts across the 
upstream rDio3 region. The single rDio3os EST along with alternatively 
spliced RT-PCR and RACE products are shown in comparison to published 
mouse Dio3os exons (Hernandez, 2005) and the 5' Dio3 genomic region 
(Dentice etal., 2007). rDio3 transcriptional (arrow) and translational (ATG) start 
sites are indicated. Locations of representative primer pairs used for 
strand-specific RT-PCR (X, Y, Z) are indicated (bold arrows), as is aT/A SNP 
used for allele-specific expression analysis. (B) Representative strand-specific 



RT-PCR of adult male rat hippocampal (X, Y) or placental (Z) total RNA using 
primers specific for opposite strand expression (reverse primer: R) or sense 
expression (forward primer: F) across the rDio3 upstream region shows 
expression is derived exclusively from the opposite strand. Controls: Positive: 
Random hexamer-primed RT-PCR (RH), which will copy potential transcripts 
from either strand; genomic DNA (DNA). Negative: Reverse transcriptase 
deficient (RT~) or template deficient (Temp - ) RT-PCR reactions; 
template-deficient (H2O) PCR reactions. 
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than an extended rDio3 5'UTR variant. We identified this poly- 
A sequence in RACE products as well. However, early reports 
of rDio3 expression based on Northern blot hybridization from 
rat brain (Tu etal., 1999) suggested the presence of alternative, 
longer (>2.1 kb) rDio3 transcripts, with later studies in human 
tissue demonstrating DI03 transcription overlapping canonical 
DI03 promoter sequences, as well as identifying a longer, 4.8 kb 
transcript that hybridized exclusively to the 5' but not 3' regions 
of DI03 and was thus suggested not to encode the DI03 protein 
(Hernandez et al., 2004). This 4.8 kb transcript would be predicted 
to overlap the 5' end of the mDio3 5'UTR-overlapping mDio3os 
splice variant (Hernandez et al., 2002) and would be expected to 
extend at least 2 kb further 5', depending on splicing profiles, 
therefore giving it further Dio3os overlap potential. However, in 
these experiments, (double-stranded) cDNA rather than allele- 
specific probes were utilized, and thus it is possible that the DI03 
promoter sequence within the extended DI03 transcript actually 
hybridized with DI03os transcripts instead. To confirm that tran- 
scription across the 5' rDio3 region identified herein in fact occurs 
exclusively in the opposite orientation to rDio3, we performed 
strand-specific RT-PCR across multiple sites, including X, Y, and 
Z (Figure 1A), in several different tissue regions/developmental 
time points (e.g., adult HP, Figure IB). As indicated in Figure IB, 
RT-PCR product was obtained only when primers complemen- 
tary to the rDio3os transcript were used. OS transcription was 
confirmed for both unspliced (e.g., across the EST BI274690, "X") 
and spliced regions ("Y"). 

IMPRINTED rDio3os EXPRESSION ARISES FROM THE SAME ALLELE(S) 
AS THE ADJACENT rDio3 TRANSCRIPT 

The rDio3os transcript had been reported to lack imprinted expres- 
sion in mouse embryos (Tierling et al., 2006) and in adult mouse 
cortex (Labialle et al., 2008). However, we have recently discovered 
that imprinting of the rDio3 transcript exhibits brain regional, 
strain, and developmental specificity (Sittig etal., 2011a), with 
the greatest differentials occurring in fetal FX, where imprinted, 
paternal expression was observed in reciprocal crosses, relaxing 
to biallelic expression in adult FX, and in adult HP, where rDio3 
expression was imprinted and derived from the maternal allele 
in one cross but biallelic in the reciprocal cross. As is custom- 
ary for antisense transcripts in imprinted loci, we expected that 
rDio3os expression would be oppositely imprinted from rDio3, 
or equally biallelic. To determine whether the rDio3os transcript 
also exhibited complex imprinted expression patterns, we identi- 
fied single-nucleotide polymorphisms (SNPs) within the rDio3os 
transcribed region between Sprague-Dawley (SD) and Brown- 
Norway (BN) rat strains. Focusing on a T/A SNP distant from 
G/C-rich or repetitive sequence regions (Figure 1A), we then 
examined expression across this SNP in animals derived from 
reciprocal crosses of these strains (S x B or B x S), where the 
first strain represents the maternal allele. We performed allele- 
specific expression analysis by strand-specific RT-PCR followed 
by direct sequencing on dissected frontal cortices and hippocampi 
from fetal and adult animals. As our prior experiments deter- 
mined that expression across the region derived exclusively from 
the OS direction, we also performed direct sequencing on prod- 
ucts isolated following random hexamer (RH) -primed RT-PCR 



reactions from these regions, as well as across several SNPs in 
placenta. 

As shown in Figure 2, rDio3os expression in brain is imprinted 
in a strain, region, and developmental-specific manner parallel- 
ing that of rDio3 (Sittig et al., 201 la). Specifically, both rDio3 and 
rDio3os show relatively biallelic expression patterns in adult FX 
from both S x B and B x S crosses, comparable to the bial- 
lelic expression reported in adult mouse cortex (Labialle etal., 
2008). In the adult HP from the B x S cross, both rDio3 (Sittig 
etal., 2011a) and rDio3os are also biallelic. In placenta, the pre- 
dominant, unspliced rDio3os transcripts, which comprise ~80% 
of total rDio3os, are also biallelically expressed (58:42 s.d. 5.6, 
n = 3), consistent with what had been observed for Dio3 in rat (Sit- 
tig etal., 2011a) and mouse (Yevtodlyenko etal., 2002), although 
others have reported preferential paternal mDio3 expression 
(Tsai etal, 2002). 

Surprisingly, however, not only is rDio3os imprinted in both 
S x B and B x S fetal FX, and in S x B adult HP as is rDio3 
(Sittig etal., 2011a), but the expression arises from the same 
allele as that of rDio3 (Figure 2). Analysis of the spliced placen- 
tal rDio3os transcripts suggest that these, also, exhibit imprinted 
expression, redolent of the conflicted status of placental murine 
Dio3 expression. The more common ("Y") variant demonstrates 
preferential paternal expression (86:14 s.d. 20, n = 3; p < 0.05 
compared to unspliced biallelic transcript), whereas a minor alter- 
native splice variant exhibits exclusively maternal expression in 
two of three placenti. In placenta, we also observe unspliced bial- 
lelic and spliced paternal expression across a downstream SNP, 
although its location distant from splice junctions and adjacent 
to a region of repetitive element sequence homology precluded 
similar examination in brain samples, which have markedly lower 
rDio3os expression levels than does placenta. Thus, unlike the 
typical sense/antisense arrangement in imprinted genetic loci, 
where the antisense transcript arises from the opposite allele 
from that of the sense transcript, in brain and potentially in 
placenta, both rDio3 and rDio3os transcripts arise from the 
same allele. 

TOTAL rDio3os EXPRESSION TRACKS rDioZ: STRAIN AND BRAIN 
REGION-SPECIFIC EFFECTS 

The standard model for sense and antisense expression in 
imprinted genomic loci is that, in general, expression of the anti- 
sense transcript occludes expression of sense transcripts, through 
direct overlap (Williamson etal, 2011) or in trans as a negative 
regulator of expression through interaction with complexes at the 
gene promoter (Shin etal, 2008). Expression of such sense and 
antisense genes is thus usually complementary, or inversely pro- 
portional. In contrast, most IncRNAs located near coding genes 
have been found to independently transcribed (Liao etal., 2011). 
Of the ~2% coregulated with coding genes, most are not found 
in the "head-to-head" orientation of the Dio3/Dio3os transcripts 
(Mercer etal, 2008; Liao etal., 2011), as currently characterized 
(Hernandez etal., 2004). As Dio3os total expression has been 
reported to be roughly correlated with Dio3 in vivo during devel- 
opment, adulthood (Hernandez etal., 2002; Labialle etal., 2008), 
and in tissue culture (Kester et al, 2006) for mouse and human, 
respectively, we wanted to determine whether this coregulation 
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FIGURE 2 | Allele-specific expression analysis demonstrates that 
rDio3os expression is imprinted in a strain, region, and developmental- 
specific manner paralleling that of rDio3. Here, various imprinted 
expression patterns are shown for fetal and adult male hippocampus (HP) 
and frontal cortex (FX), respectively. rDio3 imprinted expression is 
reproduced from our previous work (Sittig etal., 2011a). Relative rDio3os 
expression across a SNP between Sprague-Dawley (SD;T allele) and 
Brown-Norway (BN; A allele) animals was determined following (a) 



strand-specific or (b) random hexamer-primed RT-PCR, followed by (c) 
direct sequencing or (d) pyrosequencing analysis (Sittig etal., 2011a). 
For each cross, the maternal strain is given first (S x B or B x S). For 
allele-specific expression, paternal:maternal (Pat:Mat) expression ratios 
are given, along with representative preferential paternal, maternal or 
biallelic sequence traces. Strand-specific allele-specific analysis of 
rDio3os samples performed in duplicate-quadruplicate, s.d., standard 
deviation. 



extended to the regional and strain-specific total expression pat- 
tern observed for rDio3 in brain. Furthermore, we aimed to 
determine if rDio3os exhibited transcriptional response to fetal 
alcohol exposure similar to that of rDio3 (Sittig etal., 2011b), 
despite its novel imprinted, homo-allelic expression profile and 
unusual genomic organization. 

Initially, we examined total rDio3os expression from micro- 
dissected rat brain regions (Figure 3). We found that despite the 
low levels of expression of this transcript we were able to achieve 
reliable amplification using concentrated cDNA samples in a semi- 
quantitative RT-PCR assay, whereas we were unable to generate 
useable standard curves or reproducible relative expression pro- 
files for rDio3os using real-time qRT-PCR amplification cocktails. 
Previously, we had found that in the adult male S x B rat brain, 
total rDio3 mRNA expression levels are lower in the HP than in 
the FX, but these regions exhibit similar total rDio3 expression 
levels in the BxS cross (Sittig et al., 201 la). We confirmed a simi- 
lar pattern for rDio3os, wherein rDio3os expression is significantly 
reduced in the male adult HP compared with that in the FX in the 
SxB cross, but is statistically indistinct between regions in the 
BxS cross. 

We then observed that, as with native expression, total rDio3os 
response to prenatal ethanol exposure mirrors rDio3 expression. 
Figure 4A illustrates RT-PCR products from representative indi- 
vidual SxB fetal brain regions, demonstrating that expression 
across the rDio3os transcript roughly correlates with relative rDio3 
expression for that individual, as measured by real-time qRT-PCR. 
Overall, we observed elevated rDio3os expression in ethanol- 
exposed (E) vs. control (C) FX in three of six female fetal samples, 
and in three of seven adult samples, by gender (Figure 4B), 




FX HP 

rDio3 £ 5 §S 

GO CO 

FIGURE 3 | Strain and brain region-specific effects on total rDio3os 
expression. As for rDio3, semi-quantitative RT-PCR from adult male cortex 
demonstrates no significant difference between animals derived from a 
BxS cross vs. an S x B cross (p < 0.84), whereas animals from the SxB 
cross exhibit significantly decreased total rDio3os expression in the 
hippocampus vs. cortex. SxB adult hippocampi also trend toward 
decreased rDio3os expression as compared with hippocampi from the 
BxS cross (p < 0.08). For comparison, total rDio3 adult male hippocampal 
RNA expression is indicated below bars (from Sittig etal., 2011a; n = 7-8). 
Samples are normalized to rat beta-actin expression levels, n = 3-4; 
extreme outliers omitted from analysis. *p < 0.05; N.S., not significant; HP 
hippocampus; FX, frontal cortex; s.d., standard deviation. 
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FIGURE 4 | Coordinate response of rDio3 and rDio3os to prenatal 
ethanol exposure. (A) Semi-quantitative RT-PCR from fetal frontal 
cortex of individual control (C) and ethanol-treated animals (E) demonstrates 
levels of rDio3os expression tracking total rDio3 expression (relative levels 
from qRT-PCR in triplicate, normalized to beta-actin; L.S.) for both unspliced 
(X) and potentially spliced (Y) transcript-specific primer pairs. All animals are 
from an S x B cross. A placental rDio3os sample shown for comparison. 
Note that 4C is an upper outlier from controls to allow visualization of 
expression product. (B) Semi-quantitative RT-PCR from female fetal or 
adult rat S x B frontal cortex (FC) from control (C) or prenatally 



ethanol-exposed animals (E). Pink/red: female; blue: male. (C) Semi- 
quantitative RT-PCR from adult male hippocampus from control (C), pair-fed 
(PF), and prenatally ethanol-exposed animals (E) demonstrates decreased 
total rDio3os expression (region Y, shown) in ethanol-exposed animals. 
Relative rDio3 protein levels (from Sittig etal., 2011 b; n = 4-6) are given 
below each category. As in fetal brain, the direction of change is the 
same as that for rDio3 although the absolute degree of change is less. 
Samples run in duplicate, normalized to beta-actin; outliers from each three 
sample group excluded. *p < 0.05, Student's f-test. s.d., standard 
deviation. 



consistent with the overall increase in rDio3 expression observed 
in ethanol-exposed FC from these respective groups (Sittig etal, 
2011b). The exponential increase in rDio3os expression between 
E and C groups together and separately by developmental stage is 
highly significant by Chi-square analysis (p < 0.0001) although the 
absolute differences in mRNA expression levels between groups are 
not significant by two-way ANOVA. In contrast, rDio3os expres- 
sion in adult male HP was decreased following ethanol exposure 
(Figure 4C), again consistent with the pattern observed for rDio3 
(Sittig et al., 201 lb). Together, these results demonstrate that total 
rDio3os expression patterns are regulated in a similar manner to 
those of rDio3, albeit with somewhat greater variability and to a 
somewhat lesser degree than the coding transcripts. 

DISCUSSION 

The present study demonstrates the coordinate expression and 
imprinting profile of a non-overlapping IncRNA/gene pair. Specif- 
ically, we show that the rat Dio3os transcript does not obviously 
overlap the rDio3 transcript itself, that it is imprinted, and that 
it is coregulated with rDio3 both at the level of total expression 
and of imprinted expression. To our knowledge, this represents 
the first example of paired sense/OS transcripts arising from the 
same allele. The original paradigm for (long) non-coding RNAs 
(IncRNAs) was that they were usually found in imprinted gene 
clusters, overlapped at least one imprinted coding transcript (in 
antisense orientation), and were themselves reciprocally imprinted 
to the protein-coding gene (O'Neill, 2005; Peters and Robson, 
2008). Subsequent identification of roles for imprinted IncRNAs 
in the regulation of imprinted expression of non-overlapping 
genes opened our eyes to the importance of these transcripts 
as direct functional entities (Nagano etal., 2008; Pandey etal, 
2008). In turn, recently published global analyses of IncRNAs has 



highlighted many non-overlapping IncRNA/gene pairs that fre- 
quently show coordinate expression profiles (Mercer etal, 2008; 
Ponjavic etal, 2009; Liao etal, 2011), with IncRNAs implicated 
in some cases in regulation of expression of the paired gene 
(Khalil etal, 2009; Ponjavic etal., 2009). Here, we have shown 
that, much like Archaeopteryx, the rat Dio3os transcript exempli- 
fies both of these paradigms, carrying features of and serving as a 
span between the broad and generalized class of IncRNAs, and the 
more specialized category of imprinted non-coding transcripts. 

The great majority of IncRNAs that have been identified within 
imprinted gene clusters are themselves imprinted, or have tissue- 
specific imprinted isoforms (Numata etal., 2010), and thus our 
observation of region-specific imprinting of rDio3os (Figure 2), 
previously reported to exhibit biallelic expression, might not be 
considered particularly unusual. Slightly more unusual is the 
genomic organization of Dio3/Dio3os (Figure 1), which at face 
value appears to be a bidirectional, head-to-head pairing of sense 
(coding) and antisense, or OS, transcripts, with no detectable 
sequence overlap in rat (this manuscript), or human (Hernan- 
dez et al, 2004), although a small 5' overlap has been identified in 
mouse (Hernandez et al, 2002). This is in contrast to the multiple 
sense/antisense pairs of coding/non-coding RNAs with significant 
overlap found within imprinted gene clusters (Katayama etal., 
2005; O'Neill, 2005). Whereas there are also many imprinted IncR- 
NAs within these clusters that do not overlap coding transcripts 
(e.g., Zhang et al., 20 1 1 ) and which may technically be in antisense 
orientation to coding genes, including ncRNAs with functional 
roles on non-overlapping transcripts such as Air (Nagano et al., 
2008) and Kcnqlotl (Pandey etal., 2008), the distance between 
these is usually much larger than that between the identified Dio3os 
and Dio3 transcripts. There does exist a possibility of transcrip- 
tional overlap between Dio3os and Dio3, as Hernandez et al. (2004) 
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have identified a putative Dio3os promoter downstream of Dio3, 
which could facilitate chromatin reorganization across the Dio3 
locus. In addition, as neither we (rDio3os) nor others (Hernandez, 

2005) have been able to identify the absolute 5' end of the Dio3os 
transcript using RACE or other methods, there also remains a like- 
lihood that the full Dio3os transcript flanks the Dio3 gene itself, 
such that transcription proceeds across Dio3 without overlap in 
the final processed products. In either case, this organizational 
structure is highly unusual for imprinted as well as non-imprinted 
IncRNA/coding RNA gene pairs (Ponjavic et al., 2009), most espe- 
cially when the latter show evidence of coregulation (Mercer et al, 
2008; Liao et al, 201 1), as do Dio3/Dio3os (Hernandez et al, 2002; 
Figures 3 and 4). 

Most unusual is our finding that imprinted rDio3os expression 
originates from the same allele as that of rDio3. Within imprinted 
clusters, IncRNAs are almost universally transcribed from the allele 
opposite to the imprinted coding transcripts in the regulatory 
domains (O'Neill, 2005). To our knowledge, rDio3/rDio3os repre- 
sent a previously unreported and novel scenario, wherein both an 
imprinted gene and its immediately adjacent non-coding RNA are 
transcribed from the same, rather than opposite, parental alleles. 
The closest comparison to this scenario occurs within the complex 
Gnas locus, where both Nespas (overlapping Nesp55) and Gnasxl 
(one of several alternative transcripts of Gnas, along with Nesp55) 
are paternally expressed, in opposite orientation (Williamson et al, 

2006) , albeit further apart than Dio3/Dio3os (Williamson etal., 
2002). However, unlike Dio3/Dio3os, the expression of Nespas 
and Gnasxl is not concordant across tissues, although both are 
expressed in heart and brain (Pasolli etal., 2000; Wroe etal, 
2000; Plagge etal, 2004). Furthermore, the Dio3 promoter has 
not been found to exhibit differential methylation (Tsai etal., 
2002), and, as mentioned above, may not itself drive Dio3os 
expression (Hernandez etal., 2002), whereas expression of Nes- 
pas and Gnasxl is regulated by a differentially methylated region 
(DMR) between them, that controls imprinted expression across 
the locus (Williamson etal., 2006). Finally, although the imprint 
of both pairs is maintained by a distant DMR, for Gnasxl this is 
downstream of Nespas (Frohlich etal., 2010), whereas imprinted 
mDio3 expression is regulated by a DMR almost 1 Mb distant 
(Lin etal., 2003). 

This novel collection of features suggests that the regulation 
of expression of Dio3/Dio3os likely differs from the standard 
paradigm for coding/non-coding genes within imprinted loci, and 
may be better represented by that observed for coregulated cod- 
ing/lncRNA pairs. In imprinted loci, expression of the non-coding 
transcript generally forestalls expression of the coding gene on 
that allele, either by direct transcriptional (stochastic) interference, 
often leading to epigenetic changes in the chromatin of overlapped 
regulatory regions (Williamson etal, 2011), or by recruitment 
of chromatin modifying factors to (Pandey etal., 2008) or direct 
interaction with (Nagano et al, 2008) the promoters of neighbor- 
ing genes. In comparison, about 20% of non-imprinted IncRNAs 
are coregulated with their neighboring coding gene when in a 
sense/antisense orientation (Katayama etal., 2005), although cod- 
ing/lncRNA pairs in all orientations have been described that 
are coordinately, inversely, or randomly regulated (Mercer et al, 
2008). Of course, coregulation of adjacent genes may simply reflect 



a similar chromatin environment, and the unique expression pro- 
file of rDio3/rDio3os may not imply a regulatory role for either 
in generating that profile. On the other hand, IncRNA members 
of coregulated gene pairs have, in many instances, been shown to 
exhibit regulatory function on the coding gene (Katayama etal., 
2005; Mondal etal., 2010). They have also, for example, been 
identified as chromatin-associated RNAs (CARs; Mondal etal., 

2010) , that positively regulate transcription of neighboring genes 
via establishing active chromatin structures through interaction 
with the chromatin. In that study, however, Dio3os was not identi- 
fied among the CARs, although another ncRNA in the Dlkl -Dio3 
locus, Meg3, which exhibits standard opposite imprint expres- 
sion to Dio3, was. Tellingly, neither has Dio3os been pulled out 
in several other large-scale studies on IncRNAs, including those 
identified as associated with chromatin modifying factors (which 
tend to decrease expression of adjacent genes; Khalil etal, 2009), 
those identified by evolutionary sequence constraints (Ponjavic 
etal., 2009), chromatin state (Guttman etal, 2009), expression 
in brain (Mercer etal., 2008) or across transcriptome analysis 
(Katayama etal, 2005), and micro-array platforms (Liao etal., 

201 1) . In addition, Dio3os, despite being characterized in the pub- 
lished literature (Hernandez et al., 2004) and listed in IncRNAdb 
and NONCODE, databases of non-coding RNAs (Amaral etal., 
2011; Bu etal, 2012), contains no expression profiling informa- 
tion in these, nor in the Allen Brain Atlas (Lein et al., 2007) and its 
low level and tissue-specific expression may preclude it from being 
identified in such screens. 

Despite the apparently elusive nature of the Dio3/Dio3os 
partnership in large-scale screens for IncRNA functionality, the 
conservation of sequence along with coregulation of total expres- 
sion between species and maintenance of tissue-specific imprint 
patterns in wild- type and developmentally substandard conditions 
described herein, suggest that coregulation is actively maintained 
and retains the possibility of a role for Dio3os in active regulation 
of imprinted Dio3 expression. As a coregulated gene pair, Dio3os 
might enhance or, contrarily, temper expression of Dio3, both sce- 
narios having been observed for coregulated IncRNA/coding gene 
pairs (Katayama et al, 2005; Mondal et al, 2010; Vance et al, 201 1). 
Multiple potential mechanisms exist for these possibilities, similar 
to those observed for other IncRNA-based regulation. Whether 
the Dio3os transcript originates within the Dio3 promoter region 
or from the flanking, 3' end of Dio3, active Dio3os transcription 
across the locus may create an "open" chromatin structure per- 
missive for transcriptional machinery action upon the Dio3 gene. 
Alternatively, the resultant chromatin structure may limit accessi- 
bility, tempering but not occluding Dio3 transcription, as would be 
expected from an imprinted antisense transcript expressed from 
the allele opposite to that of gene expression. Direct transcrip- 
tional overlap might also regulate the post-transcriptional stability 
of Dio3 by dsRNA or siRNA formation via pre-processed Dio3os 
transcripts spanning Dio3, or an as-yet-unidentified Dio3os 5' 
exon within the Dio3 3'UTR (Hernandez, 2005). Active Dio3os 
transcription and presence of the Dio3os transcript itself adjacent 
to the Dio3 promoter, or overlapping it (Hernandez etal., 2002; 
Dentice et al., 2007), may prevent binding of an inhibitor protein, 
thus serving to activate Dio3 expression, or limit binding of a tran- 
scription factor or Pol II recruitment, thus moderating expression. 
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Another possibility is that, as for several other IncRNAs, the Dio3os 
RNA itself may be functional, homologously pairing directly with 
the Dio3 promoter region and recruiting activating or inhibitory 
factors, or regionally interacting with chromatin modification fac- 
tors to epigenetically alter the chromatin landscape, facilitating or 
tempering Dio3 expression. Any of these models offer the format of 
an additional layer of transcriptional regulatory control for Dio3, 
which, while possibly being an evolutionary relic, may in fact be 
necessary to effect proper control of total Dio3 expression levels, 
and thus of circulating thyroid hormone levels in the developing 
and adult animal, in a temporal and regional-specific manner. 

These studies demonstrate that the Dio3os transcript is present 
in rat, and that its structure is generally conserved with that in 
mouse. As in mouse, total rDio3os levels generally track those of 
rDio3. Our detailed analyses further demonstrate that total expres- 
sion similarities for the rDio3 gene and the rDio3os IncRNA are 
maintained and are similarly influenced by developmental stage, 
brain region, strain background, and prenatal insult (ethanol 
exposure). 

We further demonstrate that rDio3os expression is also 
imprinted, as is common for IncRNAs within imprinted clusters. 
As for total expression, imprinted rDio3os expression patterns vary 
in line with those of rDio3. Although this could simply be a sit- 
uation of coregulation of these genes, it is tempting to speculate 
that rDio3os expression may be positively regulating rDio3 expres- 
sion, as has been described for other non-imprinted lncRNAs/- 
adjacent genes (Katayama et al., 2005; Mercer et al., 2008; Mondal 
etal, 2010). 

Future studies including identification of rDio3os transcrip- 
tional overlap of rDio3 regulatory regions, as in mouse and 
humans, and knockdown and truncation of the Dio3os transcript 
in neuronal tissue will be required to evaluate the requirement 
of Dio3os expression in Dio3 transcriptional regulation. In com- 
parison with what is known for other IncRNA/coding gene pairs, 
the direction of Dio3o5-mediated regulation of Dio3, if any, will 
guide further studies addressing the highly unusual mono-allelic 
expression of these genes and the mechanism of their coreg- 
ulation or sequential regulation. In turn, understanding these 
mechanisms may aid in addressing malfunctions in Dio3 expres- 
sion patterns, among others, and their phenotypic consequences, 
such as are observed following prenatal ethanol exposure (Sittig 
etal, 2011b). 

MATERIALS AND METHODS 
ANIMALS 

Animal procedures and tissue collection and processing were as 
described (Sittig et al., 201 la,b) Animal procedures were approved 
by the Northwestern University Animal Care and Use Committee. 
Adult SD and BN male and female rats (70-85 days of age, Harlan, 
Indianapolis, IN, USA) were mated ( n = 6 SD females, n = 8 BN 
females) to obtain reciprocal Fl hybrid offspring, B x S and S x B, 
with first letter of the maternal strain first followed by the first letter 
of the paternal strain. For the fetal study, dams were sacrificed by 
decapitation on gestational day 21 (G21). Additional pregnant 
dams (BN n = 7; SD n = 8) were allowed to give birth and rear 
litters. One or two males and one or two females from each litter 
were sacrificed following behavioral testing (Sittig etal., 2011a,b) 



after 60 days of age. Prenatal ethanol (E) exposure and maternal 
diet procedures were performed as described previously (Sittig 
etal., 2011b). Briefly, pregnant females were assigned to a diet 
group on G8: E dams received an ethanol-containing (5% w/v, 35% 
ethanol-derived calories) liquid diet between G8 and G21 (Lieber- 
DeCarli '82; Bio-Serv, Frenchtown, NJ, USA). Pair-fed (PF) dams 
received an amount of isocaloric liquid diet (Lieber-DeCarli '82; 
Bio-Serv, Frenchtown, NJ, USA) that matched the paired E dam's 
diet consumption on the previous day. Liquid diets were replaced 
with lab chow on G21, while control (C) dams received lab chow 
and water ad libitum from Gl to G21. 

TISSUE COLLECTION 

Pregnant dams were sacrificed by decapitation on G21 between 
1000 and 1200 hours, as previously described (Sittig etal., 
20 1 1 a,b) . Fetal heads were collected directly into RN Alater reagent 
(Ambion, Austin, TX, USA) and kept at room temperature for 
24 h before brain regions were dissected and transferred into fresh 
RNA/afer. Fetal frontal cortices were dissected using a microscope 
according to the Atlas of Prenatal Rat Brain Development: A/P 
(1.0-2.2; M/L 0.0-1.5; D/V 0.0-1.0; Altaian and Bayer, 1995). 
Fetal hippocampal dissections included the entire structure. Adult 
offspring were sacrificed by decapitation 2 weeks after the final 
behavioral test. FX (A/P 5.2-2.7; M/L 0-3.3; D/V 9.0-5.0) and 
HP were immediately dissected (Paxinos and Watson, 1997). The 
left half of each structure was placed on dry ice and the right 
half was placed into RNA/afer. Tissues were stored at — 80°C 
until use. 

RNA ISOLATION AND SEMI-QUANTITATIVE RT-PCR 

Total RNA extraction was performed using Trizol reagent (Life 
Technologies, Gaithersburg, MD, USA) according to the manufac- 
turer's protocol. Genomic DNA was removed using the TURBO 
DNA-/ree kit (Applied Biosystems, Foster City, CA, USA). DNased 
RNA ( 1 |otg) was reverse transcribed using the Promega ImPromll 
Reverse Transcription kit (Promega, Madison, WI, USA) and RH 
primers. Final product was resuspended to 100 (Jtl as per manu- 
facturer's instructions. For strand-specific RT-PCR, DNased RNA 
(0.3 u,g) was reverse transcribed using rDio3 sense-strand (F) or 
OS-specific (R) primers (0.6 nM). Two microliters of final product 
(10 ul) was used directly in subsequent PCR amplification proce- 
dures, with addition of 0.2 nM of additional RT primer and 0.5 nM 
of respective Forward or Reverse primer for second strand syn- 
thesis and amplification. Controls included reactions without 
reverse transcriptase (RT - ) or template. For semi-quantitative 
analysis, RH-primed total RNA samples were amplified to obtain 
signal within the exponential phase of the PCR reaction as fol- 
lows: rDio3os: 60° C annealing temperature, 37 cycles; beta-actin: 
1:4 dilution of cDNA, 55° C annealing temperature, 29 cycles, and 
analyzed by gel electrophoresis followed by digital imaging (Kodak 
Gel Logic 200, Rochester, NY, USA) and relative densitometry 
(Adobe Photoshop 12.0.1, San Jose, CA, USA). rDio3os expres- 
sion signals were normalized to beta-actin. Primer pairs: rDio3os 
"X": F, 5'-CTTGGAGGGCCTGGCATTAAC; R, 5'-AAGACACT- 
GGCACTACTGGC; "Y": F, 5'-AACTTTCTCGACCAGAAACC- 
GC; R, 5'-TAGTATAGGAGTCCGATGGC; "Z": F, 5'-AAGCTGG- 
TTAAGGGTGGAGC; F* 5'-TTGCAACTTGAGCCCTGAGGG; 
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R, 5'-TACACCATTGCCACCACCGACTGC; beta-actin: F, 5'- 
GCTCCTCCTGAGCGCAAGTA; R, 5'-CTCCTGCTTGCTGAT 
CCACAT. 

rD/oJos CHARACTERIZATION AND SEQUENCING 

Putative rDio3os sequence was identified through BLAST 
(http://blast.ncbi.nlm.nih.gov/Blast.cgi) alignment of the EST 
BI274690 with rat chromosome 6 genomic contig reference assem- 
bly NW_047772.1 and mouse chromosome 12 reference assembly 
NC_000078.5. PCR primers for amplification of putative rDio3os 
transcript were designed manually based on mouse Dio3os exons 
in splice variant ESTs AY077459 and AY238181. To identify 5' 
and 3' transcript ends, RACE was performed using the Invitrogen 
GeneRacer kit (Life Technologies, Grand Island, NY, USA). PCR 
amplified products and rDio3os splice variants were extracted from 
agarose gels, purified and sequenced using the Children's Hospital 
of Chicago Research Center (CHCRC) Core Facility (Chicago, IL, 
USA) or the Northwestern University Feinberg School of Medicine 
(FSM) Genetics and Genomics Core Facility (Chicago, IL USA). 



For SNP identification, BN, SD, and S x B genomic DNA was 
amplified and sequenced across the region (156500-160800 bp). 
Reported analyses are using T/A (SD/BN) at 158338 bp. For allele- 
specific expression analysis, primers flanking the SNP(s) between 
SD and BN rat strains were used for amplification of cDNA derived 
from RH-primed RT-PCR or strand-specific RT-PCR (above), 
with a second round of nested PCR as necessary to generate suffi- 
cient product for sequencing. Ratio of allele-specific transcripts 
were determined by direct measurement of sequence traces in 
both forward and reverse directions, normalized to genomic DNA 
allelic ratio (50:50). Additional primers used for nested PCR and 
sequencing are available upon request. 
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